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LONG-TERM  GOALS 

•  Improve  our  understanding  of  the  ocean  basins  for  scientific  research  and  Naval  operations. 

OBJECTIVES 

•  Improve  global  marine  gravity  maps  by  a  factor  of  2  in  deep  ocean  areas  and  a  factor  of  4  in  the 
Arctic  areas  and  on  shallow  continental  margins. 

•  Use  these  improved  gravity  maps  along  with  a  global  compilation  of  soundings  to  refine  a  30- 
arcsecond  bathymetry  model  (SRTM30PLUS). 

•  Prepare  the  next  generation  of  scientists  for  ocean  research. 

APPROACH 

1 .  Modify  waveform  retracking  software  to  improve  the  range  precision  of  radar  altimeter  data  from 
CryoSat,  Envisat,  and  Jason-1. 

2.  Construct  and  assess  the  accuracy  of  new  global  marine  gravity  models  based  on  all  available  radar 
altimeter  profiles. 

3.  Construct  global  bathymetric  models  at  1 -minute  and  30-arc  seconds  using  the  new  marine  gravity 
models  and  all  publicaly-available  bathymetry  sounding  data. 

WORK  COMPLETED 

We  have  completed  the  development  of  wavefonn  retracking  software  for  CryoSat,  Envisat,  and 

Jason- 1  and  have  retracked  all  the  existing  data.  This  involved  the  development  of  a  new  analytic 

model  for  a  SAR  wavefonn  for  the  open  ocean  similar  to  the  Brown  model  waveform  used  for 

standard  radar  altimetry  [Brown,  1977].  This  research  has  been  submitted  for  publication  [ Garcia , 
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Sandwell,  and  Smith,  2013].  Using  all  these  data  along  with  the  Geosat  and  ERS-1  altimetry  data  we 
have  published  a  new  global  gravity  model  (V2 1 . 1  available  at 

ftp://topex.ucsd.edu/pub/global_grav_lmin).  This  work  was  published  in  the  August  2013  issue  of 
The  Leading  Edge  [ Sandwell  et  al,  2013].  In  June  2013,  the  Jason-1  satellite  completed  its  406-day 
geodetic  phase.  We  have  prepared  a  beta  version  of  the  global  gravity  model  (V22. 1)  that  includes  all 
the  Jason-1  altimetry  for  it  full  1 1 -year  lifetime  as  well  as  34  months  of  CryoSat-2  data  (Figure  1). 
These  results  were  presented  at  two  scientific  meetings  -  the  IAG  meeting  in  Potsdam  Gennany  and 
the  Living  Planet  Meeting  in  Edinburgh,  Scotland.  The  presentations  and  publications  are  available  at: 
ftp://topex.ucsd.edu/pub/sandwell/cryosat/Potsdam_sandwell.ppt 
ftp://topex.ucsd.edu/pub/sandwell/cryosat/LP_sandwell.ppt 
http://topex.ucsd.edu/sandwell/publications/142.pdf 
http://topex.ucsd.edu/sandwell/publications/144.pdf 

RESULTS 

Marine  gravity  and  bathymetry  are  foundational  data,  providing  basic  infrastructure  for  military, 
scientific,  economic,  educational,  and  political  work.  Naval  operations  require  accurate  gravity  models 
for  inertial  navigation  and  fire  control,  and  accurate  bathymetry  to  assess  navigational  hazards  in 
uncharted  areas.  In  addition,  coastal  bathymetry  is  used  for  improving  models  of  tides  and  currents. 

The  accuracy  of  the  global  marine  gravity  field  depends  primarily  on  two  factors  -  spatial  track  density 
and  altimeter  range  precision.  Recently  three  new  non-repeat  altimeter  data  sets  are  becoming 
available.  (1)  CryoSat  has  been  collecting  data  over  the  world’s  oceans  for  the  past  3.5  years.  While 
its  main  mission  is  to  map  the  changing  topography  of  the  ice  sheets  it  is  also  revolutionizing  our 
understanding  of  marine  gravity.  (2)  In  October  2010,  the  Envisat  radar  altimeter  was  placed  in  a  new 
partly  drifting  phase  orbit  to  conserve  fuel.  Although  the  spacecraft  failed  in  April  2012,  it  was  able  to 
collect  1.5  years  of  data  along  this  new  ground  track.  (3)  In  May  of  2012,  the  Jason- 1  radar  altimeter 
was  maneuvered  from  its  10-day  repeat  track  to  a  406  day  repeat  track  to  begin  its  new  geodetic 
mission  phase.  In  June  2013,  the  Jason-1  satellite  failed  after  1 1  years  of  data  collection.  The  failure 
occurred  just  4  days  after  the  satellite  completed  its  geodetic  phase.  This  has  resulted  in  a  high  spatial 
density  track  coverage  at  a  track  inclination  that  complements  the  high-density  data  being  collected  by 
CryoSat-2  (Figure  1). 

The  largest  error  source  in  the  construction  of  marine  gravity  models  from  satellite  altimetry  is  due  to 
errors  in  picking  the  arrival  time  of  the  return  radar  echo.  We  have  developed  a  new  approach  to 
retracking  altimeter  waveforms  by  assuming  that  the  rise-time  parameter,  which  represents  the 
roughness  of  the  ocean  surface  from  ocean  swell  (SWFI),  varies  smoothly  along  the  satellite  ground 
track.  The  results  of  the  20  Hz  range  precision  for  all  the  non-repeat  altimeters  are  provided  in  the 
Table  below.  Note  the  improved  range  precision  of  the  newer  altimeters  on  CryoSat,  Envisat  and 
Jason- 1  with  respect  to  the  older  altimeters  on  Geosat  and  ERS-1.  These  improvements  in  range 
precision  combined  with  the  increased  track  density  are  resulting  in  dramatic  improvements  in  the 
accuracy  of  the  marine  gravity  field.  A  blind  test  of  the  altimeter-derived  gravity  was  performed  by 
ConocoPhillips  using  their  proprietary  terrestrial  gravity  measurements.  The  rms  difference  between 
1.4  million  data  points  from  the  northern  Gulf  of  Mexico  is  1.60  mGal.  The  rms  difference  over  an 
area  of  semi -permanent  ice  cover  is  3.77  mGal.  Most  of  the  error  at  lower  latitude  is  due  to  errors  in 
the  E-W  component  because  CryoSat  tracks  are  mainly  oriented  N-S.  Jason- 1  provides  better  E-W 
coverage  due  to  its  lower  orbital  inclination. 


2 


Table  -  Altimeter  Noise  at  20  Hz  (mm) 


Altimeter 

2  m  SWH 

6  m  SWH 

Geosat 

57.0 

105.4 

ERS-1 

61.8 

111.8 

Envisat 

51.8 

88.6 

Jason- 1 

46.4 

64.2 

CryoSat  -  LRM 

42.7 

71.1 

CryoSat  -  SAR 

49.7 

110.9 

Standard  deviation  of  retracked  20  Hz  height  estimates  for  two  significant  wave  heights  (SWH).  The 
data  are  from  a  region  of  the  North  Atlantic  with  relatively  high  sea  state.  The  values  represent  the 
median  of  thousands  of  estimates  over  a  range  of  SWH 

Over  the  past  two  years,  we  have  completed  the  following: 

1)  Developed  a  new  analytic  model  for  a  SAR  altimeter  waveform.  The  shape  of  the  return  radar  echo 
from  an  ocean  surface  is  the  triple  convolution  of  the  source  time  function,  the  effective  area  of  the 
ocean  illuminated  versus  time,  and  the  ocean  surface  roughness  function  (i.e.  significant  wave  height). 
When  the  CryoSat  altimeter  is  operated  in  the  SAR  mode  the  form  of  the  area  of  the  ocean  illuminated 
versus  time  is  very  different  from  a  standard  altimeter  leading  to  a  sharper  return  pulse.  We  have 
developed  an  analytical  functional  fonn  for  this  return  pulse  shape  which  enables  the  use  of  standard 
least-squares  software  methods  for  rapidly  retracking  the  vast  amount  of  waveform  data  (~600  million 
waveforms/yr.) 

2)  Developed  and  refined  software  for  retracking  all  waveforms  from  Geosat,  ERS-1,  Envisat,  CiyoSat 
and  Jason.  Prior  to  this  investigation,  there  were  no  software  tools  for  two-pass  retracking  of  the 
newer  altimeter  data.  Over  the  past  two  years  we  have  developed  codes  to  read  the  raw  waveform  data 
from  each  satellite,  reformat  the  data  into  a  common  wavefonn  data  structure,  and  retrack  the  data 
using  the  2-pass  retracking  methods  [Sandwell  and  Smith,  2005]  that  are  optimized  for  gravity  field 
recovery.  Examples  of  along-track  slope  profiles  and  resulting  gravity  field  around  Hawaii  are  shown 
in  Figure  1 . 

3)  Developed  a  new  global  gravity  model  with  an  accuracy  of  1.60  mGal  in  the  Gulf  of  Mexico.  Using 
these  retracked  altimeter  profiles  from  all  available  non-repeat  orbit  satellite  altimeters,  we  have 
constructed  a  new  global  gravity  anomaly  model  (V2 1 . 1  available  at 

ftp://topex.ucsd.edu/pub/global_grav_lmin).  The  method  of  construction  is  identical  to  our  previous 
approach  [ Sandwell  and  Smith,  2009]  so  the  only  change  is  the  addition  of  the  new  data.  We  have 
assessed  the  accuracy  of  the  model  using  a  shipboard  gravity  grid  produced  by  Edcon  Co.  The  rms 
difference  between  the  satellite  gravity  and  Edcon  gravity  is  1.60  mGal  (Figure  2.).  Based  on  these 
results,  we  expect  the  gravity  accuracy  to  be  better  than  1.5  mGal  as  CryoSat -2  collects  additional  data 
and  we  refine  the  gravity  field  construction  methods. 

4)  Development  of  a  slope  correction  to  improve  gravity  accuracy  over  the  crests  of  large  seamounts. 
Recently  we  have  determined  that  the  standard  approximation  to  locate  the  reflection  point  of  the  radar 
altimeter  footprint  is  resulting  in  errors  of  several  milligals  over  areas  of  very  high  sea  surface  slope 
such  as  large  seamounts  and  trenches.  We  are  currently  assessing  the  magnitude  of  this  error  and  are 
developing  a  “slope  correction”  method  to  further  improve  the  gravity  accuracy  over  areas  having  very 
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large  anomalies  (>~  100  mGal).  This  slope  correction  will  be  implemented  over  the  third  year  of  the 
investigation. 


Figure  1.  Along  track  sea  surface  slope  from  CryoSat,  Envisat,  and  Jason-1.  The  along-track  slopes 
from  all  satellites,  both  ascending  and  descending  directions,  are  combined  using  a  biharmonic 
spline  method  to  produce  north  and  east  sea  surface  slope  maps  and  ultimately  gravity  anomaly 

(lower  right). 
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Figure  2.  (upper)  Satellite  gravity  anomalies  from  VI 8  (Geosat  and  ERS-1  only)  and  V21  (all 
altimetry  data)  plotted  versus  shipboard  gravity  from  Alaminos  Canyon,  Gulf  of  Mexico  (EDCON- 
PRJ,  Inc.).  The  accuracy  of  the  EDCON  data  is  0.5  mGal  so  discrepancy  is  due  to  errors  in  the 
altimeter-derived  gravity,  (lower)  coherence  between  satellite  gravity  and  EDCON  gravity  over  the 
area  of  the  Alaminos  Canyon.  The  coherence  falls  to  0.2  at  a  wavelength  of  20  km  for  VI 8  and  a 
wavelength  of  13  km  for  V22.  Adding  the  new  altimeter  data  from  Envisat,  CryoSat-2,  and  Jason-1 
provides  improvement  in  the  14  km  to  40  km  wavelength  band. 


IMPACT/APPLICATIONS 


Gravity  Applications: 

inertial  guidance  of  ships,  submarines,  aircraft,  and  missiles 

planning  shipboard  surveys 

mapping  seafloor  spreading  ridges  and  microplates 

continental  margin  structure 

Bathymetry  Applications: 

identification  of  navigational  hazards 

tsunami  propagation  and  inundation  models 

tide  models  and  tidal  friction 

coastal  tide  model  improvements 

ocean  circulation  models 

tidal  role  in  ocean  mixing 

understanding  seafloor  spreading  ridges 

identification  of  linear  volcanic  chains 

education  and  outreach  (i.e.  geography  of  the  ocean  basins) 

law  of  the  sea 
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TRANSITIONS 


These  global  gravity  anomaly  and  bathymetry  products  are  used  throughout  the  Navy  and  defense 
industry.  For  example,  the  National  Geospatial-Intelligence  Agency  used  our  global  marine  gravity 
grid  for  the  Earth  Gravitational  Model  2008  (EGM2008;  Pavlis,  et  al.,  2008).  The  gravity  field  is  used 
to  identify  "Red  Dots"  for  survey  planning  as  part  of  the  Submarine  Navigation  Improvement  Program 
(SNIP)  James.Braud@navy.mil.  SRTM30PLUS  is  being  used  in  the  Google  Maritime  Globe  for  the 
DoD  and  Intel  community  (paul.k.sartorius@lmco.com). 

RELATED  PROJECTS 

This  research  is  co-funded  by  ConocoPhillips  Co.  and  the  National  Science  Foundation.  In  addition, 
we  have  recently  received  funding  from  NGA  to  improve  coastal  tide  models.  All  four  sponsors  will 
benefit  from  improved  marine  gravity  and  bathymetry. 

REFERENCES 

Brown,  G.S.  (1977),  The  average  impulse  response  of  a  rough  surface  and  its  application.  IEEE 
Transactions  on  Antenna  and  Propagation,  AP-25(1):  67-74. 

Pavlis,  N.K.,  Holmes,  S.A.,  Kenyon,  S.C.,  and  Factor,  J.K.,  2012.  The  development  and  evaluation  of 
the  Earth  Gravitational  Model  2008  (EGM2008),  J.  Geophys.  Res.,  117,  B04406, 
doi:10. 1029/201 1JB008916. 

Sandwell,  D.T.  &  Smith,  W.H.F.,  (2005).  Retracking  ERS-1  Altimeter  Wavefonns  for  Optimal 
Gravity  Field  Recovery,  Geophys.  J.  hit.,  163,  79-89,  2005. 

Sandwell,  D.  T.,  and  W.  H.  F.  Smith  (2009),  Global  marine  gravity  from  retracked  Geosat  and  ERS-1 
altimetry:  Ridge  segmentation  versus  spreading  rate,  J.  Geophys.  Res.,  1 14(B0 14 1 1),  1-18. 

PUBLICATIONS 

Garcia,  E.,  D.  T.  Sandwell,  W.  H.  F.  Smith.  Retracking  CryoSat-2,  Envisat,  and  Jason- 1  Radar 
Altimetry  Waveforms  for  Improved  Gravity  Field  Recovery,  Revised:  Geophysical  Journal 
International,  June  15,  2013. 

Sandwell,  D.  T.,  E.  Garcia,  K.  Soofi,  P.  Wessel,  and  W.  H.  F.  Smith,  Towards  1  mGal  Global  Marine 
Gravity  from  CryoSat-2,  Envisat,  and  Jason- 1,  The  Leading  Edge,  August,  2013. 


6 


